Mass, radius, and composition of the outer crust of nonaccreting cold neutron stars 
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The properties and composition of the outer crust of nonaccreting cold neutron stars are studied 
by applying the model of Baym, Pethick, and Sutherland, which was extended by including higher 
order corrections of the atomic binding, screening, exchange and zero-point energy. The most recent 
experimental nuclear data from the atomic mass table of Audi, Wapstra, and Thibault from 2003 is 
used. Extrapolation to the drip line is utilized by various state-of-the-art theoretical nuclear models 
(finite range droplet, relativistic nuclear field and non-relativistic Skyrme Ifartree-Fock parameter- 
izations). The different nuclear models are compared with respect to the mass and radius of the 
outer crust for different neutron star configurations and the nuclear compositions of the outer crust. 



I. INTRODUCTION 

In 1971 Baym, Pethick, and Sutherland (BPS) deter- 
mined in their classical paper [l[ the equation of state 
and the sequence of nuclei appearing in the outer crust 
of cold nonaccreting neutron stars. They used a droplet 
model of Myers and Swiatecki Q for the masses of the 
nuclei, the only input informationen needed in this cal- 
culation. With time, more and more masses of so far un- 
known nuclei became available. Also theoretical models 
vifith increased precision were developed. Several publi- 
cations followed [3, 0, S 0] , which updated the results of 
BPS by applying the newest experimental data and theo- 
retical mass tables available at that time. In the present 
publication we take a selection of the theoretical nuclear 
mass tables and results presented in Q and analyze the 
impact of the different models on the mass radius rela- 
tion and the compositions of the outer crust for different 
realistic neutron star configurations. 

Our study of the outer crust is based on the BPS- 
model, extended by including higher order corrections. 
The only input needed are the masses of the nuclei. 
Therefore, we take the most recent experimental data 
available, the 2003 atomic mass table from Audi, Wap- 
stra, and Thibault 0- To obtain nuclear masses, 
the atomic electron binding energies were subtracted 
with an empiric formula as given in For nuclei 

with unknown masses, five different theoretical nuclear 
models were used. BSk8 is a Skyrme-Hartree-Fock- 
Bogolyubov (Skyrme- HFB) calculation [P]. SkM* is a 
Skyrme-Hartree-Fock model with a BCS pairing force 
|l(j |. FRDM stands for the most recent version of the 
finite-range droplet model fll'l . NL3 is a well known pa- 
rameterization of the relativistic mean field model [12| . 
The chiral Lagrangian used for the set Chiral is built on 
the nonlinear realization of the chiral SU(3)xSU(3) sym- 
metry as motivated from quantum chromodynamics [l3| . 
Besides the Chiral model, the nuclei were calculated al- 
lowing for axial deformations. The nuclear models used 
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are described in more detail in [5|. 

The outer crust of nonacreeting cold neutron stars is 
built of the groundstate of cold {T = 0), neutral and 
beta-equilibrated matter at densities between 10^ and 
10^^ g/car". In this density region the distance between 
atoms is small compared to their spatial extension which 
leads to complete ionization of the atoms. The nega- 
tively charged electrons separate from the positive nu- 
clei and can be described as a free Fermi-gas in first 
approximation. The nuclei arrange in a way to mini- 
mize the Coulomb interaction between themselves and 
the surrounding electron gas and thus form a solid (body 
centered cubic) lattice. 

The total pressure of the system P has a contribu- 
tion from the electrons Pe and a negative pressure com- 
ing from the lattice: P = Pe + \WLnM^ where Wl is 
the lattice energy and the number density of nuclei. 
The total mass of the nuclei Wn enters in the baryon 
chemical potential = {Wn + + ^Me)M, with the 
mass number A and the charge number Z . The electron 
chemical potential /.tg is fixed by the condition of charge 
neutrality. Thus the baryon chemical potential at given 
pressure is only a function of the charge and mass num- 
ber. To find the groundstate, one has to look for the 
nucleus (A, Z) with minimal baryon chemical potential. 
A more detailed description of the BPS model is given in 
the original paper [l| , or can be found in |5|] . 

So far, the electrons were considered as free particles. 
In the present calculation, three higher order corrections 
to this approximation are included |14| . First the screen- 
ing or Thomas-Fermi energy, which describes deviations 
of the electron distribution from uniformity due to the 
negatively charged ions. Secondly, the exchange energy, 
to take into account the Fermionic nature of the electron- 
electron interaction. Furthermore the zero-point motion 
of the nuclei in the lattice was incorporated, also to verify 
the stability of the solid crystal. The latter correction to 
the energy density turns out to be even smaller than the 
first two and is negligible with respect to the equation of 
state or the composition of the outer crust. 

With the lowest energy per nucleon of all elements, 
^^Fe represents the ground state nucleus up to densities of 
10^ g/cTP? . At higher densities the electrons become rel- 
ativistic and begin to change the position of equilibrium. 
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Electron capture takes place, leading to more and more 
neutron rich nuclei. Exotic nuclei with large neutron ex- 
cess appear, stabilized against /3-decay by the filled elec- 
tron Fermi sea. Finally, the baryon chemical potential 
exceeds the neutron mass. Now the nuclei are so neu- 
tron rich that some neutrons arc not bound any more 
and start to drip out of the nucleus. The neutron drip is 
reached and the inner crust begins. 

The equilibrium nucleus is determined by an interplay 
of the strong, short ranged nuclear interactions present 
in the nuclei and the electromagnetic forces between the 
charged particles. In this regard the outer crust is an 
fascinating astrophysical site, where theoretical nuclear 
models determine macroscopic, in principle observable, 
properties of neutron stars. 

The internal structure of a neutron star is described by 
the Tolman-Oppenheimer-Volkoff (TOV) equation p^ : 
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with the radius r, the mass m, the gravitational con- 
stant G, and the speed of light c. With the knowledge 
of the equation of state, i.e. the energy density e{P) as 
a function of the Pressure P, the mass and the radius 
corresponding to one certain central pressure Pq can be 
calculated. Usually one starts the numerical integration 
of the TOV-equation in the center of the star, i.e. r = 0, 
m = 0, and P = Pq. But the TOV-equation also allow 
to begin at an arbitrary point inside the neutron star, 
r — Rq, m — Mo, and P — Pq, as far as realistic values 
for Rq, Mq and Pq are considered. This allows to calcu- 
late the mass and radius of the outer crust for a chosen 
inner neutron star configuration {Rq, Mq and Pq), with- 
out any assumptions about the equation of state of the 
inner region below the outer crust. 

The mass and the radius are the most elementary and 
important measurable quantities of a neutron star. Pre- 
cise mass and radius measurements can rule out equa- 
tions of state which are in conflict with the experimental 
data, as most recently demonstrated in the case of the 
X-ray binary EXO 0748 - 676 [H. Astrophysical obser- 
vations are a source of knowledge about the properties of 
dense hadronic matter. In this context it is interesting to 
study the consequences of the theoretical nuclear models 
on the mass and radius of the outer crust, which is the 
topic of this investigation. 

On the other hand, the composition of the outer crust 
is important for heat transport and for Ohmic dissipation 
of the magnetic field, which will be explored in the last 
section. Furthermore, the composition and the equation 
of state are necessary to determine the eigenmodes of 
acoustic oscillations or toroidal shear modes [isj . 



II. MASS AND RADIUS OF THE OUTER 
CRUST 

Table |T] lists the mass and the radius calculated for 
three different equations of states, which are based on the 
theoretical mass model BSk8 and SkM*. Additionally 
one equation of state for BSkS is considered, where effects 
from the lattice (and the higher order corrections to it) 
were neglected. The initial pressure for the beginning of 
the outer crust is chosen to be Pq = 6.6 x 10^^ dyne/cm^, 
a value equivalent to the neutron drip of these equations 
of state. 



TABLE I: The mass AM and thickness AT? of the outer crust 
of a nonaccreting cold neutron star with radius Rq = 12 km 
and mass Mo = 1.4 Mq of the inner region (inner crust to- 
gether with core), for three different equations of state. The 
initial pressure for the beginning of the outer crust is fixed to 
Po = 6.6 X 10^^ dyne/cm^ 



EoS 


AM [10-^ Mq] 


AR [km] 


BSkS 


3.090 


0.4509 


SkM* 


3.088 


0.4408 


BSkS without lattice 


3.093 


0.4666 



Compared to the total radius of a neutron star of 10 
to 20 km the deviations in the radial extensions of the 
outer crusts are found to be quite small. In this respect 
the theoretical nuclear models only play a minor role in 
the mass radius relation. However it is surprising, that 
the differences of the two nuclear mass tables BSkS and 
SklVI* are almost as important as the inclusion of the 
lattice. 

Next, the dependence of the mass and thickness of the 
outer crust on the radius Rq and mass Mq of the inner 
region will be analyzed. The equation of state which 
belongs to the mass table FRDIVI is being used. Figure 
[1] depicts the results for initial masses of Mq — 1.0 — 2.0 
Mq and initial radii of i?o = 10 — 20 km. 

The mass as well as the thickness of the outer crust in- 
creases with Rq for all initial masses. On the other hand 
neutron stars with large Mq have smaller and lighter 
outer crusts. For the considered values of Mq and Rq the 
outer crust mass AM ranges from 7.8 x 10~® to 5.6 x lO^'', 
the radius AR from 0.14 to 2.57 km. For Mq = 1.0 Mq 
the outer crust spans between 4.7 and 11.4 % of the total 
radius R + AR, for Mq = 2.0 Mq between 1.4 and 4.7 
%. In summary, the outer crust is the lighter and smaller 
the more compact the neutron star is. 



III. COMPOSITION 

In the following the composition of the outer crust 
for different neutron star configurations and equations 
of state is examined. The distribution of the fraction 
of nuclei on the total number of nuclei Ntot for given 
mass number A and charge number Z is shown in Figure 
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H The theoretical nuclear models NL3, FRDM and the 
more exotic model Chiral were selected for comparison. 

The sequence of groundstate nuclei with increasing 
density for FRDM and NL3 follows mainly along the 
magic neutron numbers TV = 50 and iV = 82 [Si], which 
leads to the bimodal distribution visible in Figure [21 In 
all three models there is only a vanishingly small amount 
of ^^Fe present, as this nucleus only appears at the low- 
est densities. On the contrary, elements which appear 
at high densities can give a significant contribution to 
the overall composition, like it is the case for ^^®Pd in 
the Chiral model. Although the groundstate nucleus are 
identical for all the three models up to 3x 10^° g/cm'^, this 
causes the visible differences in the total compositions. 
These results clarify that differences in the sequence are 
likely to be amplified in the composition. 

With NL3 a mean mass number (A) = 105.7 and mean 
charge number (Z) = 34.5 is found for the Mq = 1.4 M0 
and i?o — 10 km neutron star, {A) — 108.3 and {Z) — 
36.0 for FRDM, and (A) = 104.8 and {Z) = 34.0 for the 
set Chiral. These values are relatively similar, in average 
the groundstate nuclei are Selenium {Z = 34) or Krypton 
(Z = 36) nuclei, with mean mass number of 105 to 108. 

The fractions in Figure [2] were calculated for two dif- 
ferent neutron stars, for Mq — 1.4 Mq and Rq = 10 km 
and a less compact neutron star with Mo = 1.0 M0 and 
Rq — 20 km. Although the crusts of these two differ- 
ent neutron stars are very different in size (see Figure [T]) 
the distributions of the nuclei are almost unaltered. This 
means that all existing nonaccreting cold neutron stars 
should exhibit a very similar composition of their outer 
crusts. Of course this conclusion is only valid, if the neu- 
tron star has actually reached the absolute ground state 
during its evolution and if one neglects contaminations 
of any fall-back supernova material. 

IV. SUMMARY 

For given neutron star configurations, the mass and 
thickness of the outer crust of nonaccreting cold neutron 
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FIG. 1: The mass AM and thickness AR of the outer crust in dependence of the inner radius Ro for different inner masses 
Mo, calculated for the mass table FRDM. The initial pressure of the outer crust is Po = 7.9 x 10^^ dyne/cm^. 
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FIG. 2: Distribution of the fraction of nuclei on the total number of nuclei in the outer crust Ntot in dependence on the mass 
number A (left side) and the charge number Z (right side) for the equations of state NL3, FRDM, and Chiral. Black bars show 
the results for a neutron star with Mo = 1.4 Mq, i?o = 10 km, gray ones for Mo = 1.0 M0, Ro = 20 km. The corresponding 
crusts are 288 m and 2,754 m thick with 1.61 x 10"^ M© and 5.59 x 10"* M©, respectively. Po = 7.9 x 10^^ dyne/cm^. 



